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摘 要 : 基于 多 源 .多 期 遥感 影像 与 DEM 提取 的 冰 湖 数据 ,分 析 了 念 


青 唐 古 拉 山 西 段 1991 一 2017 年 冰 湖 不 同 维度 


的 分 布 变 化 ,在 此 基础 上 研究 了 冰 湖 数量 、 


看 积 变化 与 气温 .降水 变化 的 关系 。 


湖 数 量 与 面积 均 呈 上 升 趋势 , 冰 湖 数量 增加 161 个 


,面积 增加 1.185 km’, 


结果 表明 :(1) 近年 来 ,研究 区 的 冰 
其 中 ,1991 一 2006 年 平稳 增长 ,2006 一 


2017 年 大 幅 增长 ,.(2) 冰 湖 面积 与 数量 在 中 高 海拔 范围 内 均 呈 增长 趋势 。 新 增 冰 湖 主 要 分 布 在 海拔 $400~5700 m, 
冰 湖 数量 占 整体 数量 的 72.41%。(3) 分 布 在 不 同 朝向 冰 湖 的 增长 速率 不 同 。 西 南 向 冰 湖 增长 较为 迅速 , 朝 北向 冰 


湖 增 长 相对 缓慢 
湖 对 气候 的 响应 更 灵敏 ,变化 更 显著 ,面积 大 的 冰 湖 7 
过 探究 念 青 
关键 词 : 念 青 唐 古 拉 山西 段 ; wee; 冰 湖 ; 气候 变化 


冰川 湖泊 是 在 冰川 退缩 消融 后 为 了 填充 侵蚀 、 
消 蚀 形成 的 洼地 而 形成 的 蓄 水 体系 ,简称 冰 湖 。 大 
多 数 冰 湖 位 于 水 文 基准 面 之 上 " ,其 生长 受 季 节 和 
年 份 等 时 间 尺 度 的 影响 ,是 气候 变化 与 冰川 退缩 共 
同 作用 所 形成 的 产物 ,因此 多 数 冰 湖 是 短暂 的 , 冰 
湖 被 认为 是 全 球 变化 的 天 然 指 示 器 之 一 ”。 在 全 球 
范围 气候 变 暖 趋势 下 ,高 亚洲 地 区 冰川 后 退 现象 显 
著 汪 ,而 冰川 后 退 常 伴随 着 冰川 湖泊 的 变化 : 。 

冰 湖 的 形态 及 其 变化 情况 不 仅 能 侧面 反映 出 
全 球 气 候 的 整体 变化 趋势 ,也 能 为 评价 冰 湖 演 决 灾 
害 做 出 预警 "。 青 藏 高 原 包含 许多 大 型 冰 磺 基地 湖 
泊 , 这 些 湖泊 近 几 十 年 来 不 断 扩 大 "* ,因此 ,青藏 高 
原 最 近 冰 湖 演 决 (Glacial Lake Outburst Flood, 
GLOF ) 的 频率 和 幅度 Pn 有 效 获取 
冰 湖 信息 并 对 冰 湖 的 状态 进行 监测 分 析 , 可 以 帮助 
ee eet 页 
测 、 预 报 。 

近年 来 国内 外 主要 对 喜马拉雅 山地 区 *”” KR 
山地 区 外、 别 珍 套 山 等 地 区 的 冰 湖 变化 进行 研 
究 讨论 ,关于 气候 与 冰 湖 耦合 研究 方面 主要 是 针对 
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Heme: 


(4) 气温 上 升 引 起 的 冰川 消融 是 冰 湖 生长 的 间接 因素 ,降水 是 冰 湖 生长 的 重要 因素 。 
变化 幅度 小 ， 
甫 古 拉 山 西 段 冰 湖 变化 对 冰 湖 演 决 灾害 的 预测 具有 重要 的 指示 意义 。 


面积 小 的 冰 
表明 冰 湖 的 动态 变化 与 气候 变化 紧密 相关 。 通 


喜马拉雅 山 ” 珠穆朗玛 峰 "”、 科 西河 流域 及 阿 
尔 泰山 "等 大 尺度 地 区 进行 研究 ,而 针对 小 尺度 地 
区 ,如 念 青 唐古拉 山地 区 的 研究 极 少 ,而 念 青 唐 古 
拉 山 作为 高 亚洲 地 区 冰 湖 主要 分 布 的 区 域 ,对 其 冰 
湖 特征 及 其 变化 进行 分 析 具 有 实际 意义 。 


1 研究 区 概况 与 方法 


1.1 研究 区 概况 

念 青 唐古拉 山 位 于 西藏 纳木错 地 区 与 林芝 地 
区 的 交界 处 , 呈 东 西 走向 , 自 西 起 四 底 斯 山脉 ,向 东 
北 延 伸 , 至 那曲 附近 与 横断 山脉 相连 ,东北 向 紧邻 
雅鲁藏布江 最 大 的 支流 拉萨 河流 域 。 海 拔 范围 在 
3000~7000 m ,平均 海拔 5500 m。 核 心 研究 区 位 于 念 
青 唐 古 拉 山 西 段 (30"34'~30?34'N , 89°51’~90°49’E) , 
研究 区 DEM 影 像 的 轮廓 范围 如 图 1。 念 青 唐 十 拉 山 
脉 西 段 的 冰川 属于 亚 大 陆 性 冰川 ,而 山脉 东 段 却 受 
印度 洋 西南 季风 影响 显著 ,降水 多 ,雪线 海拔 低 , 冰 
川 分 布 集中 , 占 整个 山岳 冰川 总 量 的 80% 以 上 ,是 
中 国 最 大 的 大 陆 性 冰川 集中 地 区 ,也 是 地 球 上 中 低 
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图 1 研究 区 位 置 示意 图 
Fig. 1 Location of the study area 


纬 地 区 最 强 的 冰川 作用 中 心 之 一 。 

念 青 唐古拉 山西 段 属于 半 干 旱 大 陆 性 气候 ,年 
均 气 温 常 位 于 0 以 下 ,年 降水 量 在 300~400 mm。 
夏季 , 当地 气候 受 温 暖 潮湿 的 季风 气 团 控制 ,主要 
从 南部 穿 过 雅鲁藏布江 山谷 。 季 风气 团 携带 大 量 
水 分 ,使 东部 念 青 唐 古 拉 山 成 为 整个 青藏 高 原平 均 
降水 率 最 高 的 地 区 ,6 一 9 月 降水 量 可 达 全 年 降水 量 
的 80% 。 夏 季 日 最 高 气温 可 达 12 % ,冰川 更 易 产 生 
剧烈 的 消融 促进 冰 湖 的 生长 ;冬季 ,西风 主导 了 环 
流 模 式 , 主 要 输送 来 自 青 藏 高 原 的 冷 干 气 团 ,导致 
冬季 降水 量 少 且 气温 低 , 最 低 气 温 为 -13.7%C。 
1.2 数据 来 源 

本 研究 所 采用 的 数据 资料 包括 遥感 影像 数据 、 
DEM 数据 和 气象 数据 。Landsat 系列 数据 来 源 于 美 
地 质 调查 局 (USGS, https://glovis.usgs.gov/app ) 和 
地 理 空间 数据 云 (http://www.gscloud.cn/) , Google 
Earth 遥感 影像 来 源 于 BIGEMAP (http://www.bige- 
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map.com/)。 在 选取 遥感 影像 时 , 需 注意 影像 成 像 时 
间 与 云 量 覆 盖 率 等 因素 的 影响 ,选取 9 一 10 月 冰 湖 
形态 稳定 的 秋季 影像 ,同时 要 挑选 云 量 较 小 的 数据 ， 
降低 云 量 对 研究 过 程 中 提取 冰 湖 信息 的 干扰 。 除 此 
之 外 ,采用 Google Earth 15 级 数据 ,分 辨 率 为 4m 的 
和 遥感 影像 作为 参考 , 评 佑 信息 提取 结果 的 精度 ”。 
遥感 影像 数据 基本 信息 见 表 1。SRTM DEM 数据 来 
源 于 美国 国家 航空 航天 局 (NASA , https://urs.earth- 
data.nasa.gov ) ,分 辩 率 为 90 mx90 m, 数 据 重 采样 后 
得 到 30 m 分辨 率 的 DEM 数据 ,基于 该 数据 派生 的 
高 程 . 坡 向 坡度 等 信息 进行 冰 湖 分 布 变化 分 析 。 
气象 数据 来 源 于 中 国 科学 院 资源 环境 数据 云 平台 
(http://www.resdc.cn/Default.aspx ) ,选取 1990—2015 
年 全 国 逐 年 平均 气温 年 降水 量 空间 插值 数据 集 。 
与 气象 站 点 无 法 获得 偏远 地 区 数据 有 旦 检测 到 的 数 
据 呈 分 散 点 相 比 ,该 数据 集 是 经 过 插值 处 理 后 得 到 
空间 平均 温度 ,获得 的 气象 数据 精度 更 高 ,可 信 度 
更 强 。 
1.3 研究 方法 

国内 外 运用 遥感 技术 提取 水 体 的 方法 较 多 , FE 
要 包括 目 视 解 译 手工 数字 化 ”基于 光谱 信息 的 效 
值 法 和 NDWI 指 数 法 ”、 基 于 机 器 学 习 的 卷 积 神经 
网 络 .决策 树 和 支持 向 量 机 等 方法 ,以 及 基于 图 像 
分 割 方法 的 “全 域 - 局 部 ”法 玻 和 面向 对 象 法 等 ”。 
除 此 之 外 ,还 有 “全 域 -局 部 ”与 NDWI 相 结合 、 
Sentinel-1 图 像 与 Landsat 8 图 像 结 合 ”” 及 半自动 提 
取 ” 等 方法 。 因 本 文 研究 区 域 范围 较 小 ,人 研究 对 
象 又 多 分 布 在 冰川 末端 , 冰 湖 面积 大 多 又 小 于 0.01 
km ,采用 自动 或 半自动 方法 提取 的 冰 湖 轮廓 并 不 
准确 ,后 期 还 需 人 工 核 对 , 故 选取 能 够 直接 获取 高 


表 1 研究 区 遥感 数据 


Tab.1 Remote sensing data of the glacial lakes in the study area 


传感器 成 像 时 间 轨道 参数 平均 云 量 /% a Ga 用 途 
Landsat 4/5 TM 1991-09-15 138/039 7 30/120 冰 湖 提取 
Landsat 4/5 TM 1996-10-13 138/039 4 30/120 冰 湖 提取 
Landsat 7 ETM+ 2001-11-04 138/039 3 15/30 冰 湖 提取 
Landsat 5 TM 2011-10-07 138/039 4 30/120 冰 湖 提取 
Landsat 7 ETM+ 2010-09-26 138/039 38.12 15/30 参考 数据 
Landsat 7 ETM+ 2012-10-17 138/039 30.3 15/30 参考 数据 
Landsat 8 OLI 2017-10-07 138/039 5.58 15/30 冰 湖 提取 
Landsat 8 OLI 2017-10-23 138/039 2.67 15/30 参考 数据 
Google Earth 2017-10-09 一 4 参考 数据 
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精度 冰 湖 信息 的 手工 数字 化 方法 。 借 助 地 理 信息 
系统 软件 完成 冰 湖 的 矢量 化 , 因 山 体 阴影 与 冰 湖 不 
易 区 分 , 故 在 提取 冰 湖 前 根据 影像 数据 的 太阳 高 度 
及 太阳 方位 角 进 行 山体 阴影 演 染 ,矢量 化 过 程 中 保 
证 误差 不 超过 一 个 像 元 值 ,将 提取 结果 与 DEM 派生 
数据 坡度 图 进行 琶 加 分 析 , 坡 度 小 于 5 的 水 体 认 定 
为 冰 湖 ,大 于 5° 的 则 被 认为 是 山体 阴影 等 干扰 项 ， 
去 除 干扰 项 后 保留 其 余 冰 湖 提 取信 息 '"。 人 参考 
Google Earth 15 级 (分 辩 率 为 4m) 的 影像 .相近 时 期 
影像 数据 和 地 形 图 辅助 判别 。 对 冰 湖 提取 结果 进 
行 空间 分 析 以 供 后 期 分 析 冰 湖 变 化 ,获取 冰 湖 的 各 
项 属性 ,包括 冰 湖 编号 、 冰 湖面 积 、 冰 湖 平均 高 程 、 
冰 湖 坡度 、 坡 向 、 冰 湖 变 化 类 型 和 冰 湖 规模 。 

假设 冰 湖 面积 的 提取 误差 符合 高 斯 分 布 ”, 则 
可 通过 冰 湖 的 周 长 确 定 像 元 个 数 .影像 的 最 大 分 辩 
率 和 随机 权重 误差 来 计算 影像 解 译 中 所 产生 的 面 
HRZ ,公式 如 下 : 


wo -学 

式 中 : u, 为 冰 湖 面积 误差 (km); p 为 冰 湖 周 长 
(m); A 为 像 元 最 大 分 辩 率 (30 m); o (0.6872) ABE 
机 权重 误差 ,假设 69% 的 像 元 会 出 现 误差 。 

影像 质量 、 像 元 分 辩 率 .几何 配 准 技术 以 及 人 
工 经 验 等 均 会 对 冰 湖 信息 提取 的 精度 产生 影响 。 
本 文 主要 侧重 探讨 念 青 唐古拉 山区 域 冰 湖 整 体 在 
1991 一 2017 年 的 变化 规律 ,而 不 是 以 单个 冰 湖 作为 
对 象 研究 逐个 像 元 的 变化 ,所 以 造成 的 误差 都 是 系 
统 性 的 。 另 外 ,在 时 间 变 化 分 析 的 过 程 中 ,边界 误 
差 会 在 6 期 数据 的 对 比 中 相互 抵消 一 部 分 ,因此 , 像 
元 配 准 的 误差 对 念 青 唐古拉 山区 域 冰 湖 整体 变化 
规律 认识 的 影响 有 限 。 

采用 参数 统计 的 线性 相关 系数 法 分 析 了 气温 、 
降水 对 冰 湖 面积 及 数量 变化 的 影响 。 线 性 相关 系 
数 又 称 为 简单 相关 系数 、Person (皮尔 森 ) 相 关系 数 
等 , 常 以 符号 > 表示 样本 相关 系数 , p 表示 总 体 相关 
系数 。 总 体 相关 系数 的 定义 公式 如 下 : 

Cov(X, Y) 
式 中 : Cov(X, Y) 是 随机 变量 X.Y 的 协 方差 ;Var@) 
和 Var( 六 分 别 代 表 对 和 YY 的 方差 。 总 体 相关 系数 是 
反映 两 变量 之 间 线 性 相关 程度 的 一 种 特征 值 ,表现 
HTE KX AX’ AX” Y’ Ye Y,D 


(1) 


Pxy =Corr(X, Y) = (2) 


别 是 总 体 和 了 一 组 简单 随机 样本 ,样本 相关 系数 的 


n 


YX, — XY, ~ Y) 
p= =) (3) 


Sa T HL, T iy 


样本 相关 系数 是 根据 样本 观测 值 计算 的 ,抽取 
的 样本 不 同 ,其 具体 的 数值 也 会 有 所 差异 。 的 绝 
对 值 越 大 表明 两 变量 间 的 相关 性 越 大 。 


2 结果 与 分 析 


2.1 总 体 分 布 特征 

通过 上 述 方法 得 到 1991 年 .1996 年 .2001 年 、 
2006 年 .2011 年 和 2017 年 的 冰 湖 信息 (图 2)。 对 研 
究 结果 进行 统计 分 析 , 得 出 1991 一 2017 年 冰 湖 数量 
和 面积 变化 特征 ( 表 2)。 结 果 表 明 ,从 1991 一 2017 
年 的 27 a 间 ,该 区 域 冰 湖 数量 增加 了 161 个 ,面积 增 
加 1.185 km ,数量 增长 率 73.51% ,但 面积 增长 率 仅 
为 22.37% ,数据 表明 这 期 间 冰 湖 数量 剧 增 ,但 冰 湖 
面积 并 没有 同步 大 幅 增 长 。 

从 冰 湖 变化 趋势 的 整体 情况 看 , 冰 湖 的 数量 和 
面积 均 呈 增长 趋势 。 其 中 ,1991 一 1996 年 冰 湖 面积 
增长 最 快 ,增长 速度 为 0.075 km-a”; 2001—2006 4 
冰 湖 面积 增长 最 慢 , 为 0.021 km?+a'; 2011—2017 4 
冰 湖 面积 增长 适中 ,但 冰 湖 数量 增长 最 快 ,增长 率 
达 18.01% ,2001 一 2006 年 冰 湖 数量 增长 与 面积 增长 
速度 一 致 ,增长 率 仅 6.59% ,为 5 个 时 段 增长 幅度 最 
小 的 时 期 。 

2.2 不 同时 期 冰 湖 变化 特征 

1991—2017 年 念 青 唐古拉 山西 段 冰 湖 存在 新 
JM 消失、 退缩 以 及 扩张 现象 ,其 空间 分 布 特征 见 图 
3。 冰 湖 变化 主要 出 现在 研究 区 的 东北 部 ,以 扩张 
冰 湖 和 新 增 冰 湖 为 主 ,分 别 占 31%、38%, 消 失 冰 湖 
数量 最 少 , 仅 占 6%。 

念 青 唐古拉 山西 段 1991 一 2017 年 新 增 冰 湖 137 
个 ,多 数 集中 分 布 于 研究 区 5500~5800 m 的 中 高 山 
区 ,新 增 冰 湖 的 面积 为 1.373 km2 ,在 冰 湖 4 类 变化 类 
型 中 其 数量 与 面积 变化 均 是 最 大 的 变化 类 型 。 新 
增 冰 湖 平均 面积 为 0.008 km? ,其 中 ,冰川 终 矿 湖 是 
新 增 冰 湖 的 主要 类 型 ,也 是 念 青 唐 古 拉 山 西 段 冰 湖 
规模 变 大 的 主要 类 型 。 图 4a 为 两 个 典型 新 增 冰 湖 
的 空间 分 布 , 冰 湖 位 置 分 别 位 于 海拔 5623 m 和 
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图 2 冰 湖 提取 结 


Fig.2 The extraction results of glacial lakes 


表 2 1991 一 2017 年 冰 湖 数量 与 面积 变化 
Tab.2 Glacial lakes changes during 1991-2017 


年 份 面积 km 数量 /个 较 前 一 时 期 面积 变化 km 较 前 一 时 期 增加 的 冰 湖 数量 /个 
1991 5.297 219 

1996 5.674 253 0.377 34 

2001 5.913 273 0.239 20 

2006 6.022 291 0.109 18 

2011 6.276 322 0.254 31 

2017 6.482 380 0.206 58 


图 3 1991 一 2017 年 念 青 唐古拉 山西 段 冰 湖 不 同 变化 类 型 
空间 分 布 
Fig.3 Distribution of different change types of glacier lakes 
of the study area during 1991-2017 


5560 m Mh , 均 是 因 冰 川 退缩 形成 的 冰川 终 磺 湖 呈 。 
念 青 唐 古 拉 山西 段 1991 一 2017 年 共 消 失 23 个 冰 
湖 ,消失 面积 为 0.328 km。 消 失 冰 湖 类 型 主要 为 冰 


磺 湖 2 汪 ,图 4b 为 两 个 典型 消失 冰 湖 的 空间 分 布 ,分 别 
位 于 海拔 为 5350 m #15380 m 处。 从 逐 感 影像 判断 ， 
消失 冰 湖 可 能 是 由 于 终 矿 堤 溃 决 ,从 而 使 冰 湖 下 干 。 

念 青 唐古拉 山西 段 1991 一 2017 年 有 90 个 冰 湖 
退缩 ,退缩 面积 为 2.902 km 。 退 缩 冰 湖 多 数 分 布 在 
海拔 5400~5600 m, DAYKA BE ZEW ANE. BI 5a 为 
位 于 海拔 5$683 m 人 处 典型 退缩 冰 湖 的 空间 分 布 , 该 冰 
湖 主 要 是 由 于 冰川 末端 的 持续 退缩 , 冰 湖 缺少 足够 
的 冰川 融 水 补给 而 缩小 呈 。 念 青 唐 古 拉 山西 段 
1991—2017 年 扩张 冰 湖 共 110 个 ,扩张 后 的 冰 湖 面 
积 为 2.119 km ,主要 分 布 在 5500~5900 m 的 中 高 山 
区 ,变化 类 型 主要 为 冰 矿 阻 奢 湖 中 的 侧 矿 阻塞 湖 与 
ATi SAFE”. Sb 为 典型 扩张 型 冰 湖 ,该 湖 位 于 
海拔 5670 m 人 处 ,属于 受 冰 川 侧 奔 蕉 阻挡 积 水 而 形成 
FAA se PAC 
2.3 不 同 海拔 冰 湖 变化 特征 

念 青 唐古拉 山地 区 冰 湖 分 布 在 海拔 4900~6100 
m 范 围 之 内 ,从 整体 看 , 冰 湖 数量 趋 于 单 峰 分 布 模式 
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图 5 典型 退缩 P IKIRE KI 
Fig. 5 The shrinkage glacier lakes and expansion of glacier lakes in the study area 


(图 6) ,峰值 主要 出 现在 海拔 5$400~5700 m ,该 海拔 
区 间 内 的 冰 湖 数量 占 整个 研究 区 的 72.41%, 冰 湖面 
积 占 总 面积 的 62.6%。 以 每 100 m 为 一 个 高 度 带 ,将 
整个 海拔 范围 分 为 11 个 区 间 进 行 分 析 , 得 到 5500~ 
5600 m 和 5600~5700 m 为 冰 湖 的 主要 分 布 区 域 , 冰 
湖 数量 与 冰 湖 面积 占 比 均 高 于 其 他 海拔 范围 。 此 
外 , 念 青 唐古拉 山 冰川 末端 海拔 趋 于 正 态 分 布 , 且 
集中 分 布 在 5400~5700 me , 这 与 该 区 海拔 5400~ 
5700 mm 的 冰 湖 变化 较为 显著 相 一 致 , 故 中 高 海拔 区 


响 也 较 小 , 冰 湖 形态 不 易 发 生变 化 ;海拔 5300~5400 
m 的 冰 湖 面积 从 1991—2001 年 呈 增 长 状态 ,2001 一 
2017 年 又 逐渐 减少 ,面积 持续 减少 至 0.582 km’, IK 
少 幅 度 为 21.8% ; 海拔 5400~5500 m、5500~5600 m 
和 5600~5700 m 的 冰 湖 在 2011—2017 年 趋 于 稳定 
增长 状态 ,增长 面积 分 别 以 7.2% .6.12% 和 4.53% 的 
平均 速率 增长 , 冰 湖 的 数量 也 逐年 增加 ,增幅 分 别 
为 24%、87.72% 和 77.05%, 因 此 ,1991 一 2017 年 念 青 
唐古拉 山西 段 区 域 的 冰 湖 在 5400~5500 m、5500~ 


域 的 冰川 消融 退缩 是 造成 念 青 唐古拉 山区 域 冰 湖 
形成 与 演变 的 主要 原因 。 
从 1991 一 2017 年 不 同 海拔 冰 湖 数量 和 面积 来 


5600 m 中 高 海拔 处 的 变化 ,是 冰 湖 面积 与 数量 逐年 
增加 的 主要 区 域 ;海拔 $800~5$900m 的 冰 湖 在 2001 
年 后 数量 及 面积 也 急速 增长 ,1991 一 2017 年 冰 湖 数 


看 ,位 于 5800~5900 m 的 冰 湖 面积 与 数量 变化 最 大 ， 
数量 增长 的 同时 面积 也 随 之 增长 。 相 对 于 中 高 海 
拔 区 域 冰 湖 的 显著 变化 ,海拔 较 低 区 域 冰 湖 变化 则 
较 小 ,海拔 小 于 5300 m 的 冰 湖 面积 变化 相对 稳定 ， 
说 明 低 海拔 区 域 的 冰 湖 距离 冰川 较 远 , 受 冰川 的 影 


量 增幅 为 73.52% , 冰 湖 面积 以 0.025 km”…(5a) 的 速 
度 增 长 ,这 表明 高 海拔 冰 湖 对 气候 变化 的 反应 更 灵 
敏 ,也 表明 高 海拔 冰川 的 退缩 为 冰 湖 的 形成 提供 了 
物质 来 源 淖 。 除 此 之 外 , 念 青 唐古拉 山西 段 冰 川 末 
端 海拔 为 5500~7000 m, 且 海拔 高 于 5700 m 的 冰川 
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6 1991 一 2017 年 念 青 唐古拉 山西 段 不 同 海拔 冰 湖 面积 与 数量 变化 


Fig.6 The number and area variation of different-elevation glacier lakes in western Nyaingentanglha range during 1991-2017 


RARE Si EK EK IME 5500~5800 m 冰 湖 
变化 较 显 著 的 结果 一 致 , 即 冰 川 退缩 是 冰 湖 数量 及 
面积 增加 的 直接 原因 。 
2.4 不 同 朝向 冰 湖 变化 特征 

念 青 唐古拉 山西 段 冰 湖面 积 在 不 同 朝 向 上 具 
有 差异 ,主要 分 布 在 南 、 东 、 东 南 和 西南 向 (图 7a)， 
面积 占 比 分 别 为 22.68% 、21.67% 、18.01% 和 
17.75%, 共 占 冰 湖 总 面积 的 80.2% ;西北 向 最 少 , 仅 
占 总 面积 的 7.2% ,其 次 为 朝 北 向 , 占 总 面积 
9.99%, 


(a) 面积 变化 北 mm 2017 年 
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1991—2017 年 念 青 唐古拉 山西 段 冰 湖 在 不 同 
朝向 上 的 数量 与 面积 变化 ( 表 3) , 冰 湖 数量 在 南 和 
西南 向 的 变化 最 为 显著 ,其 中 朝 南 向 变化 最 大 的 时 
期 为 1991 一 1996 年 ,数量 增加 了 近 一 倍 (图 7b); 西 
南 向 变化 最 大 的 时 期 也 在 1991 一 1996 年 ,数量 增加 
了 近 2 倍 , 北 和 西北 向 的 冰 湖 变化 幅度 较 小 ,尤其 是 
朝 北 癌 ,基本 上 趋 于 稳定 ,这 与 该 地 区 冰川 变化 一 
致 , 念 青 唐古拉 山西 段 冰 川 在 北向 上 退缩 程度 弱 
于 南 向 ,冰川 融 水 对 其 附近 的 冰 湖 面积 及 数量 扩 
aK OTR m 
人 北 = 2011—20174F 
eit tm 20012008 


m 1996—20014F 
== 1991—19964F 


东 


图 7 1991 一 2017 年 念 青 唐 古 拉 山 西 段 不 同 朝向 冰 湖 分 布 


Fig.7 The area variation of different-aspect glacier lakes in western Nyainqentanglha range during 1991-2017 
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#3 1991 一 2017 年 不 同 朝向 冰 湖 数量 与 面积 变化 
Tab.3 Variation of different-aspect glacier lakes changes during 1991-2017 
1991 年 1996 4E 2001 年 2006 4E 2011 年 2017 4 

iil TR gwm WR mwn UA et TR 数量 个 TR gw Pate 
dt 0.001 1 0.005 3 0 0 0.005 4 0.005 3 0.009 6 
东北 0.47 19 0.0203 7 0.047 13 0.062 16 0.04 12 0.056 13 
东 1.232 37 1.48 19 1.496 27 1.358 36 1.113 24 1.601 27 
东南 1.203 55 0.604 44 0.846 54 1.268 63 0.928 48 0.442 55 
南 1.928 69 2.556 105 2.903 117 2.861 108 2.721 123 3.007 157 
西南 0.386 18 0.885 51 0.539 45 0.361 38 1.243 68 1.111 68 
西 0.055 12 0.09 17 0.038 9 0.08 16 0.172 35 0.178 35 
西北 0.021 8 0.035 了 0.044 8 0.026 10 0.053 9 0.079 19 


2.5 不 同 规模 冰 湖 变化 
湖泊 是 陆地 表面 具有 一 定 规模 的 天 然 注 地 的 
蓄 水 体系 ,是 湖 盆 .湖水 以 及 水 中 物质 组 合 而 成 的 


由 表 4 可 知 ,1991 一 2017 年 念 青 唐古拉 山西 段 
冰 湖 面积 增加 主要 为 <0.01 km 与 0.01~0.05 km? 的 
冰 湖 。 面 积 <0.01 km 的 冰 湖 数量 最 多 ,日 生长 速度 


自然 综合 体 ” ,通常 分 布 在 坡度 较 小 的 洼地 或 深 
谷 , 地 形 起 伏 复杂 的 高 原 山 地 地 貌 , 以 冰 缘 地 貌 为 
代表 ,难以 形成 中 大 规模 湖泊 。 而 冰川 湖泊 就 分 布 
在 典型 的 冰 缘 地 貌 地 带 , 是 冰川 作用 后 在 平坦 洼地 
形成 的 自然 蓄 水 体 ,规模 往往 较 小 。 在 分 析 过 程 中 
将 冰 湖 面积 分 为 5 种 规模 ， RD aa 


湖 数量 及 面积 情况 ( 表 4)。 念 青 唐古拉 山西 段 最 小 
冰 湖 面积 低 至 0.0002 km’, 最 大 冰 湖 面积 可 达 0.332 


km2, 面 积 平均 值 为 0.017 km。 其 中 , 绝 大 多 数 冰 湖 
为 面积 小 于 0.05 km 的 小 型 冰 湖 ,面积 小 于 0.01 km? 
的 冰 湖 数量 占 比 超过 冰 湖 总 数量 的 50%。 

分 析 发 现 ,面积 <0.01 km 的 冰 湖 与 0.01~0.05 


最 快 ,数量 增长 速度 为 27.6 个 (5a)', 该 规模 冰 湖 总 
面积 除 在 2006 年 有 所 减少 ,其 余 时 期 均 呈 增长 状 
态 ; 面 积 为 0.01~0.05 km 的 冰 湖 数量 逐年 递增 , 且 
面积 随 着 冰 湖 数量 增加 稳定 增长 ,增幅 为 0.633 
km ,五 种 不 同 规模 冰 湖 中 增幅 最 大 。 以 上 两 种 规 
模 冰 湖 是 念 青 唐古拉 山西 段 冰 湖 面积 增加 与 数量 
增多 的 主要 因素 ,说 明 该 区 域内 小 规模 的 冰 湖 演化 
更 为 激烈 。 面 积 为 0.05~0.1 km 的 冰 湖 面积 1991 一 
2001 年 均 呈 增长 状态 ,但 2001 一 2017 年 有 增 有 减 ， 
面积 峰值 出 现在 2011 年 ,为 1.347 km。1991 年 起 
面积 为 0.1~0.2 km 冰 湖 数量 与 面积 均 呈 下 降 趋 势 
直至 2011 年 ,2011 年 后 面积 为 0.1~0.2 km? 冰 湖 数 


km 的 冰 湖 数量 占 冰 湖 总 数量 的 85% 以 上 ;面积 在 
0.1~0.2 km 与 >0.2 km 冰 湖 数量 仅 占 5% ,总 的 来 看 
冰 湖 面积 愈 大 数量 愈 少 。 从 不 同 规模 面积 总 和 来 
看 ,面积 在 0.01~0.05 km: 冰 湖面 积 总 和 占 比 最 大 ， 

可 达 30% 以 上 ,其 次 为 面积 为 0.10~0.20 km? HKG, 


表 4 1991- 2017 年 念 


量 及 面积 呈 增 长 状态 ,但 未 回 到 1991 年 的 峰值 ; 面 
F>0.2 km 的 冰 湖 较为 稳定 ,面积 增幅 仅 为 0.087 


km ,数量 增幅 为 0, 且 该 规模 冰 湖 均 位 于 低 海拔 ,说 
明 大 面积 、 低 海拔 的 冰 湖 形态 较为 稳定 , 受 冰川 退 
缩 及 气候 因素 的 影响 小 ,有 一 定 维稳 能 力 。 即 小 型 


青 唐古拉 山西 段 不 同 面积 规模 冰 湖 变化 


Tab. 4 The number and area variation of different-sized glacier lakes in western Nyaingentanglha range during 1991-2017 


1991 年 1996 4 2001 年 2006 4E 2011 年 2017 4E 
TTS 数量 个 “面积 数量 个 面积 数量 个 面积 数量 个 面积 。 数量 个 面积 。 数量 /个 ”面积 
/km /km /km /km /km /km 
<0.01 117 0.532 144 0.693 161 0.814 183 0.723 209 0.901 255 0.981 
0.01~0.05 75 1.611 81 1.79 82 1.8334 79 1.881 81 1.857 97 2.244 
0.05~0.10 13 0.864 16 1.127 19 1.318 15 1.042 20 1.347 15 1.037 
0.10~0.2 12 1.792 8 1.125 8 1.251 10 1.381 10 1.555 11 1.636 
>0.2 2 0.497 4 0.939 3 0.696 4 0.994 2 0.615 2 0.584 
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冰 湖 的 面积 波动 大 ,对 气候 的 反应 更 敏感 。 相 反 ， 
大 型 冰 湖 对 气候 变化 的 反应 迟钝 ,保持 原状 的 能 
更 强 不 易 产 生变 化 。 
2.6 冰 湖 变化 特征 与 气象 因素 影响 分 析 

众多 指示 气候 的 因子 中 气温 和 降水 能 直接 反 
映 气 候 变 化 , 且 是 影响 冰 湖 变化 的 主要 因子 ,降水 
量 与 气温 的 变化 会 影响 冰川 的 生长 与 退缩 ,冰川 退 
缩 消融 与 降水 都 是 形成 冰 湖 的 物质 来 源 。 

从 念 青 唐古拉 山西 段 年 均 气 温 变化 图 可 以 看 
出 (图 8a) ,研究 区 气温 变化 波动 较 大 ,但 年 均 温 度 
均 低 于 0% 。 其 中 ,年 均 最 低温 度 -3.38 °C ih SE 
1997 年 ,最 高 温度 -1.07 % 出 现在 2010 年 。 虽 然 气 
温 波 动 幅 度 大 ,无 固定 变化 规律 ,但 该 区 域 年 均 温 
度 趋势 线 1990 一 2015 年 总 体 呈 持续 上 升 趋势 ,这 与 
施 雅 风 等 人 研究 得 出 小 冰期 盛世 (17 世纪 ) 以 来 中 
国 西部 山区 气温 增加 的 结论 一 致 ” 。 近 年 来 , 念 青 
唐古拉 山西 段 气 温 逐 渐 上 升 导致 了 冰川 的 剧烈 退 
A) ,冰川 的 消融 为 冰 湖 的 形成 提供 了 物质 来 源 ， 


y= -0.0032x?+0.1269x-3.174 
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EAS SCE FE AT Be 5 TAT AS HM 2 RA BX 

从 念 青 唐古拉 山西 段 逐 年 年 降水 量 数据 来 看 
(图 8b) ,年 降水 量 均 高 于 400 mm ,气温 变化 幅度 比 
降水 变化 幅度 大 ,年 降水 量 集中 在 400 ~ 800 mm, 
整体 与 气温 变化 趋势 一 致 ,年 降水 量 也 呈 上 升 趋 
势 , 与 气温 变化 不 同 , 降 水 量 的 增幅 小 ,上 升幅 度 较 
平缓 。 冰 湖 的 增加 对 降水 量 的 变化 反应 敏感 ,新 生 
冰 湖 和 冰 湖 的 扩张 可 能 跟 降水 有 关 , 这 与 本 文 研究 
不 同 冰 湖 类 型 中 ,新 生 冰 湖 和 扩张 冰 湖 数量 最 多 的 
结果 一 致 。 

运用 皮尔 森 相 关系 数 法 分 析 气 温 降水 与 冰 湖 
变化 之 间 的 相关 性 ,由 相关 性 分 析 结 果 可 知 ( 表 5)， 
气温 与 冰 湖 面积 、 数 量 的 相关 系数 的 绝对 值 都 很 
小 ,相关 性 并 不 显著 ,气温 对 冰 湖 变化 的 影响 是 间 
接 的 ,在 其 他 气候 条 件 不 变 的 情况 下 ,气温 升 高 会 
导致 冰川 融化 进而 促进 冰 湖 面积 的 增长 。 降 水 与 
冰 湖 面积 数量 之 间 的 相关 系数 绝对 值 接近 0.6, 相 
比 与 气温 因素 ,降水 因素 与 冰 湖 变化 之 间 的 相关 性 


Y= -0.4643x2+15.153x+552.44 
R=0.1 

一 一 年 降水 
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图 8 1990 一 2014 年 念 青 唐 古 拉 山 西 段 年 均 气 温和 年 降水 量变 化 


Fig. 8 Variations of annual mean temperature and annual precipitation from 1990 to 2014 


表 5 气温 ,降水 与 冰 湖 面积 ,数量 的 相关 性 分 析 
Tab. 5 Correlation between temperature, precipitation 


and glacier lakes changes 


气温 降水 面积 数量 
气温 1 
降水 -0.491 1 
面积 0.033 0.584 1 
数量 -0.163 0.591 0.969" 1 


注 :*#* 表 示 在 0.01 水 平 ( 双 侧 ) 上 显著 相关 。 


更 为 显著 ,说 明 降水 是 影响 冰 湖 变化 的 重要 因素 。 


3 结论 


<i 


(1) 1991 一 2017 年 念 青 唐 古 拉 山西 段 的 冰 湖 持 
续 增长 ,1991 一 2006 年 平稳 增长 ,面积 增加 至 6.022 
km’, 2006 一 2017 年 涨幅 增 大 ,面积 增加 至 6.482 
km?; 27 a 来 冰 湖 数量 增加 161 个 ,面积 增加 1.185 
km’, 

(2) 冰 湖 面积 与 数量 在 中 高 海拔 均 呈 增长 趋 
势 ,5300~5400 mm 的 低 海拔 冰 湖 面积 呈 减 小 趋势 。 
新 增 冰 湖 主要 分 布 在 海拔 5400~5700 m, 冰 湖 数量 
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占 整体 数量 的 72.41%, 说 明 高 海拔 处 的 冰 湖 对 气候 ciers, Glacial Lakes and Glacial Lake Outburst Floods. Monitoring 
变化 的 响应 更 为 敏锐 ， 同时 受 冰 川 退缩 的 影响 更 大 。 and Early Warning Systems in the Hindu Kush- Himalayan Re- 

一 | bey ` . = ion: Bhutan[M]. Kathmanddu, Nelal: International Center for Inte- 
(3) 不 同 朝向 的 冰 湖 分 布 与 增长 速率 不 同 。 西 ie 
i mi vk 湖 增 长 = $ 讯 束 胡 北 i 冰 湖 增 KE M K grated Mountain Development, 2001. 
0 y ET H y ht ZF 小 有 Ss 
STs eee is BNE [10] Ageta Y, Iwata S, Yabuki H, et al. Expansion of glacier lakes in re- 
= = =I Ay = EY Vb» Aas 
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Spatial changes of the glacial lakes in the western Nyaingentanglha Range 


TAO Jing', ZHAO Wenji, WANG Xu’, XUE Jin’ 
(1. College of Resources Environment and Tourism, Capital Normal University, Beijing 100048, China; 
2. School of Geography and Information Engineering, China University of Geosciences, Wuhan 430074, Hubei, 


China; 3. Shenzhen Longhua District Development Research Institute, Shenzhen 518110, Guangdong, China) 


Abstract: Based on the multi-source remote sensing image and DEM data, we analyzed the distribution and analysis 
of different dimensions of glacial lakes in western Nyainqentanglha from 1991 to 2017. Then we analyzed the 
relationship between changes in ice lakes and temperature and precipitation. The research results showed that the 
number and area of ice lakes in the study area had increased recently. They increased steadily from 1991 to 2006 
and expanded from 2006 to 2107. The number of ice lakes increased by 161, and the area by 1.185 km in 26 
years. The area and number of ice lakes are increasing in the middle-and high-altitude range. The newly-added 
ice lakes are mainly 5400-5700 m above sea level, and the number of ice lakes accounts for 72.41% of the total 
number. The distribution and growth rates of glacial lakes with different aspects are different. The glacial lake on 
the southwest slope grew rapidly, while the north slope’ s glacial lake grew slowly. Glacier melting caused by rising 
temperatures is an indirect factor for ice lake growth, and precipitation is an essential factor for ice lake growth. 
Ice lakes with smaller area grades have a more sensitive response to climate and more significant changes. Ice 
lakes with larger area grades have stronger stability and smaller changes. The dynamic change of ice lakes is 
closely related to climate change, and it is instructive for predicting glaciers change and ice lake collapse 
disasters in western Nyaingentanglha. 


Keywords: western Nyaingentanglha Range; remote sensing; glacial lake; climate change 


